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We discuss dense states of QCD matter formed in high-energy hadronic and heavy-ion collisions 
from the point of view of statistical physics of non-equilibrium processes. For this sake, we first 
propose a formulation of the dynamical entropy of dense QCD states in the "saturation regime" 
leading to a color glass condensate (CGC). The statistical physics description amounts to describe 
the modification of the color correlation length with energy as a compression process for which 
nonequilibrium thermodynamic properties are applicable. We derive an expression of the dynamical 
entropy in terms of the rapidity evolution of the unintegrated gluon distributions in the colliding 
nuclei, verifying suitable positivity and irreversibility properties. We extend this approach to the 
' initial pre-equilibrium (glasma) state of an heavy-ion collision. It allows for a definition of the 

Q , initial entropy before the evolution towards the hydrodynamic regime as a function of the glasma 

correlation length and an overlap parameter characterizing the low-momentum spectrum of the 
glasma state. This initial entropy, by extension to the TV = 4 SYM theory, is then matched as the 
key input parameter to the strong coupling evaluation of thermalization towards the hydrodynamic 
regime based on the AdS/CFT correspondence. It thus allows to cast a bridge between the weak 
and strong coupling phases of an heavy-ion reaction. 

S 1 

1 '« I. INTRODUCTION 

J> 1 It is commonly believed that statistical physics concepts may be relevant to describe the development and outcome 
t-H , of ultra-high energy collisions between incident hadrons or nuclei. For instance, the multi-production of particles 
during such a high-energy reaction has been related to the entropy produced by the collision. A statistical physics 
framework seems especially relevant for the central distribution of multiplicities in ultra-relativistic heavy-ion collisions, 
where relativistic hydrodynamics, assumed during a transient evolution of the quark-gluon plasma (QGP), give a nice 
account [l[ of the observable features of the reaction. 

However, the theoretical difficulty remains in relating these features to the basics of Quantum Chromodynamics 
(QCD), being the field-theoretical framework where, ultimately at least, those statistical physics features, and in 
particular the entropy release, could be properly defined. Indeed, one faces both difficult statistical-physics and field- 
theoretic problems, since the initial reaction medium is both far from equilibrium and submitted to a time-dependent 
QCD evolution in a strong-coupling regime. 

Recent interesting steps forward made in this direction are the following. In the context of strong gauge coupling 
studies using the AdS/CFT correspondence [2j, it was shown [3j that starting from various initial conditions, the 
thermalization leading to a hydrodynamic behavior of a strongly-coupled gauge theory plasma essentially depends 
on an "initial entropy" factor defined within the gauge/gravity duality framework using involved notions of black 
hole (BH) thermodynamics. The hydrodynamic regime is identified with the long-time behavior of a receding BH 
In this case, the entropy release during the evolution towards the hydrodynamic regime may be associated with 
the difference between the final hydrodynamic entropy and the abovementionned "initial entropy" . This encouraging 
result has yet some limitations. It applies to a maximally supersymmetric theory in the strong coupling regime, while 
one is in practice dealing with QCD with weak-coupling initial states in the context of the Color Glass Condensate 
Q and the formation of an initial glasma [§]. Hence the problem is open to make the matching between a weak 
coupling initial phase and the strong coupling scenario of thermalization. Our goal is to know how to define properly 
the "initial entropy" for dense QCD states already at weak coupling. 

In this context, a recent paper [7] introduces the notion of a thermodynamical entropy associated with the production 
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of gluons in the saturation regime of dense CGC initial states in proton- proton collisions 1 . It is motivated by the 
relation between the saturation scale (denoted by Q s in transverse momentum) with a temperature of thermally 
produced hadrons and a dynamically generated mass, itself also related to the saturation scale Q s . This, together 
with the expected relation between entropy and the inclusive gluon production in a dense-dilute hadron-hadron 
collision, leads to an estimate of the entropy of a dense CGC initial state in the framework of the Golec-Biernat 
Wiisthoff model [|| for the QCD unintegrated gluon distribution (UGD). 

Ref.[7] is a stimulating approach, since it proposes a macroscopic approach to the entropy of a CGC ensemble of 
gluons, in a thermodynamic context. It may be hindered by the assumptions of thermodynamic equilibrium, such 
as the existence of a temperature and an entropy which are not yet well defined in principle in a perturbative QCD 
framework which seems to govern far-from-equilibrium processes where these notions are not standard. However, the 
notion of dense QCD states, which is basic in the CGC approach, has specific features which may lead to a convenient 
statistical physics description of the high-energy collision processes without making assumptions on thermodynamic 
equilibrium. In some sense it would be useful to have a microscopic definition of entropy, analogous to the Boltzmann 
statistical mechanic approach, based on the QCD dynamics of the CGC. We will propose such a notion of a dynamical 
entropy for dense QCD states in the present paper. 

In the next section [Til we w iU give the definition of a "dynamical entropy" for dense QCD states. Then in section 
Hill we will apply it to models of the CGC medium, and check the matching with the macroscopic approach of Ref. [7| ■ 
In the next section ITVl we will propose an extension of the dynamical entropy approach to the initial glasma phase of 
heavy- ion collisions and analyze its matching with the thermalization scenario based on the AdS/CFT correspondence. 
A final section [V] is devoted to conclusions with both a summary of results and outlook. 



II. DYNAMICAL ENTROPY OF CGC STATES 



Dense QCD states may be present in various physical situations. In this section, we will consider the CGC medium. 
It plays a role in Deep Inelastic Scattering (DIS) in the "small-x" range where it represents the quantum states of the 
target as seen by the virtual photon when the Bjorken variable x ~ Q 2 je Y <C 1, Q 2 being the photon virtuality and 
Y is the total rapidity range of the final states. It also appears in the description of "dense-dilute" collisions such as 
proton-nucleus at very high energies. Later on, we will extend our analysis to the initial medium in ultra-relativistic 
heavy-ion reactions obtained by the collision of two boosted nuclei in their center-of-mass frame. 

These dense QCD states may be theoretically described within QCD at weak coupling, through the nonlinear energy 
evolution starting from a dilute partonic state. The QCD evolution with total rapidity Y increases the density of 
partons and finally reaches the saturation regime leading to a CGC state. It is characterized among other features 
|5(, by a limiting transverse size R S (Y) ~ 1/Q S (Y), where Q S (Y) is the rapidity-dependent saturation momentum. 
R s (Y) is in fact the color correlation length at given rapidity. 

As energy increases, the parton density becomes high enough such that individual parton branching is compensated 
by their recombination. This leads to many interesting properties which can be evaluated still in the (resumed at 
leading or next- leading logarithms) weak coupling regime of QCD. Among them "geometric scaling" states [l(| that 
the unintegrated gluon distribution <f>(k; Y) of a CGC medium essentially depends only on one variable implying the 
saturation scale. It verifies a scaling property as a function of a single variable u = k 2 R 2 (Y), within some controlled 
approximation [111 E3 > namely 

4>{k,Y) d 2 k ~ 4>{u = k 2 R 2 ) R 2 S d 2 k . (1) 

The balance between parton branching and recombination at saturation and geometric scaling will be the main 
property allowing for a statistical mechanic approach, despite the probable absence of thermal equilibrium. Let us 
consider thi s Q CD compensation mechanism between branching and recombination from the statistical physics point- 
of-view on |13l420l | (choosing some useful references inside an abundant recent statistical physics literature, see e.g. 

The CGC distribution of gluons can be obtained by the evolution of a partonic state (essentially gluons) from a 
given rapidity Y\ to a dense medium at higher rapidity Y%, characterized by a spectrum in transverse momentum k 
defined by the UGD (f)(k, Y^). Using the formulation of classical 2 statistical physics, the partons forming the CGC 
medium may be interpreted as a stationary state, since, at each value of rapidity gluon branching and recombination 



1 A different approach to entropy production in heavy-ion collisions can be found in Refs.Q. 

2 Though obtained through a quantum field theory framework, the system of partons generated by the QCD evolution has been proven 
to be described at leading logarithm approximation as a classical branching and recombination process. 
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compensate each other in the saturation regime. The evolution in rapidity Y\ — ¥ Yi inside the dense regime can be 
considered as an (a "priori nonequilibrium) process driven by a "dynamical parameter", being the energy increase 
exp Yi/2— »exp Y2/2 of the partonic system. 

Indeed, in the statistical physics framework, a similar situation is met when a distribution of stationary states is 
defined through a probability distribution p^ t&t (z\ A), where z is a general notation for the phase space, see [ill H3] 
(and for a brief reminder, Appendix |B|). By changing the dynamical parameter Ai — > A2, the system evolves in time 
leading to a final distribution which can be compared to a new stationary distribution 'P Stat (z; A2). A key property 
of this nonequilibrium process is the Hatano-Sasa identity involving both •p Stat (, 2 :; \ 12 ), see Eq. (|Bl[) . 

Indeed, we are to show that the branching-recombining mechanism, responsible for the geometric- scaling prop- 
erty of the unintegrated gluon distribution (f>(k,R s ) ~ <p(kR s (Y)) gives rise to a transverse-momentum probability 
distribution verifying the analogous of the Hatano-Sasa identity (|B1I) . 

We are led to the following identification (see Appendix [Bjl : 

z~*k, A -» Y, P stat (z; A) dz ~* V(k, Y) <fk, (2) 
where, using ([!]). the transverse momentum probability distribution of gluons in the CGC is defined as 

the last relation coming from the geometric scaling form (ffj. This definition ensures the normalization condition 

P(k; Y) d 2 k = J P(u) du = 1 . (4) 
It is now straightforward to write the identity 

~cW {losV{k ' Y)} dY ) = I v{k ' Y2)d2kxcxp '[ log {wM} 

J / y 2 

= J P(k,Y x ) d 2 k = l, (5) 

where the left-hand side correlator (• • -) Y is defined by averaging over the probability distribution in the final state 
at Y2 and making use of the probability normalization relation (j4]) at Y± . 

As we shall demonstrate further on, Eq. ([5]) appears as the QCD version of the Hatano-Sasa identity (|B1|) for a CGC 
ensemble of states. It is a general property of QCD saturation and the CGC [|| that the gluons organize themselves in 
cells of typical saturation scale size. This is made quite explicitly, when introducing the geometric scaling formulation 
of (JT|) . Hence, the rapidity increase is responsible of the shrinkage of the color correlation length which appears as a 
compression process understood as the mean-field result of the denser gluonic medium. It gives rise to a modification 
of the probability distribution of gluon momentum which verifies the relation ([5]). This is explicit from (fT]), where the 
compression is manifest in the variation of the saturation scale R s (Yi) — > R S (Y2) < R s (Yi). 

Thanks to this interpretation, and following the proposed identification ^ , we introduce the the notion of dynamical 
entropy as follows. We define the dynamical entropy density of a CGC parton medium at rapidity Y2, coming from 
the QCD evolution in rapidity Y±— tY^ inside a transverse area of order of the initial color correlation size R 2 s (Yi) by 

We note that in the statistical physics literature formulas similar to (|6|) appear also (see e.g. as the relative 

entropy of a probability distribution w.r.t. a reference one (here V(k, I2) compared to V(k,Yi)). It amounts to 
quantify the "amount of disorder" provoked by the rapidity evolution Y\ — > Y% of the CGC medium. This notion 
is well adapted to the CGC medium of gluons, for which the competition between branching and recombination is 
expected to increase with the densification of gluons. 
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Let us list the role and properties of the dynamical entropy formula ([6]). 

1) Positivity. By construction, the identity ([5]), using the well-known [2l| Jensen inequality < (e x ) , for 

7 (fc,y 2 ) 

>(fc,Yi) • 

n^i)/y- 2 - \ 



any probabilistic distribution here identified to X = — log ^[^'y^ , leads to the positivity condition 



= ( log Zn'^i ) - 1< « ( <-xp { l< 



= 0, (7) 

Y 2 



_P(fc,Yi) 

for all Y 2 > 3^i- Hence, the familiar positivity condition is obtained for any increase in total rapidity. 

ii) Geometric scaling case in terms of non- equilibrium thermodynamics. Considering now the geometric scal- 
ing property ([I}, the identity ([5]) may also been discussed in thermodynamic terms, following the discussion and 
properties of Jarzynski identity, see Appendix 1X1 and (|A1I) . Indeed, using the geometric scaling version of ([1}, the 
identity ([5]) can be rewritten in geometric-scaling variables where for short notation i? 12 = Ri(Y 12 ). one gets 

where we reintroduced the saturation momentum scales Qg(li,2) = (Y\ o). 

The relation (jSJ) can be put in tight analogy with the Jarzynski identity [13] of the thermodynamics of far-from- 
equilibrium processes (see Appendix [X| which writes 

e -W/T\ =e -(F -F A )/T (Q) 
/ B 

where the average is made over different realizations of a process driving rapidly an equilibrium state A at tempera- 
ture T to an out-of-equilibrium state B, which thermalizes towards a new equilibrium state C (keeping the driving 
dynamical parameter constant at its new value) at the same temperature. 

The Jarzinsky identity relates the stochastic distribution of thcrmodynamical works in the process A —> B to the 
free energy balance AF between the two equilibrium states A — y C. Interestingly, the amount of dissipative work 
yyDiff = w „ AF during the process A^B is then related to the entropy production AS* = ((W) - AF)/T > 0, if 
the state B is able to relax towards the temperature T, keeping the driving parameter constant. 

Comparing now ([8]) to (HJ), one realizes that the expression logR^/Ri corresponds to the logarithm of a ratio of the 
available phase-space for dimension R\ vs. R\. It is the free energy change of the particle of an ideal gas contained inside 
a 2-dimensional "box" when the size changes by compression R±— tR® < Ri- Hence, a thermodynamic interpretation 
of the QCD relation (JSJ is that the modification of the total rapidity Y\ — > Y2 induces a modification of the CGC 
ensemble of states with reduced saturation size R2 (or equivalently increased mean momentum Q2), resulting in an 
entropy eventually generated by further on relaxation. One is led to the analogy 

B V{Rlk 2 ) T 
, R\ AF 

fog — k 

6 R\ T 

, T(k,Y 2 ) W-AF W Diff 
log ; 



V(k,Yt) T T 

k,Y!) 

V{k,Y 2 ) 



^SMH) -as'. cud 



Y 2 



where AS is, in the thermodynamic context, the entropy production (by gluonic degree of freedom) due to the 
compression R\ R2 when the system relaxes to a state at the same temperature T as the initial one but within the 
restricted domain size R 2 . 

The last line of (|10l) is suggestive of a special relation between the dynamical entropy of a CGC state S and 
the entropy production after relaxation AS. We shall confirm this analogy for Gaussian gluon distributions, by an 
explicit calculation in the next section. 



iii) Dynamical entropy of a physical CGC state. The dynamical entropy formula ([6]) remains written for a 
unit gluon degree of freedom within the color correlation area and for an arbitrary rapidity evolution Y\ — > Y 2 within 
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the QCD saturation domain. In order to go to CGC states one deals with in high-energy collisions between hadrons 
and/or nuclei, one is due to extend the formalism to those physical cases. 

Using the thermodynamic analogy, the right-hand side of the identity ([8} can be interpreted as coming from the 
standard free energy of one gluon d.o.f. inside of an ideal gas of gluons confined in a radial transverse box of saturation 
size R s , as expressed in the second raw of (|10[) . In order to determine the overall dynamical entropy that we may 
associate to a physical CGC initial state (e.g. hadronic or nuclear one), we will assume that the probability distribution 
of reference V(k, Yi) considers an initial "transverse cell" at rapidity Y\ where R\ is a typical color correlation length 
for hadronic matter (proton or nucleus) at rest. In fact, this reference distribution will not play a role for large enough 
rapidity. 

Then, we consider the boosted CGC medium when Y% — >Y%, R% — tR^, where Y 2 is the rapidity at which one observes 
[e.g. by deep-inelastic or dense-dilute hadronic scattering) the corresponding CGC state. In this framework, by 
addition of the individual contributions ©, the total dynamical entropy density dSr/dy for a boosted target CGC 



state will be the sum over all degrees of freedom, namely color multiplicity N% — 1, standard [22j gluon occupation 
number ~ 1/ ATrN c a s in longitudinal coordinate space. For the transverse degrees of freedom, one has to take into 
account the average number R 2 ^ / R\ of "transverse cells" at initial rapidity Y\ , where Rt is the hadronic or nucleus 
target size, and finally the average number fi of gluonic degrees of freedom inside a cell which will be determined later 
on using the Gaussian models for the UGD's. 
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dS _ N 2 - 1 R T yRl ^ R2 

Gluon correlation effects beyond the ideal gas approximation may induce correction factors, which we do not discuss 
in the present study. 

III. APPLICATION TO GAUSSIAN CGC MODELS 

We want now to apply our concept of dynamical entropy cf. ^ to simple CGC models in order to describe its main 
properties in a concrete way. Another motivation is to confront our definition of dynamical entropy for a physical CGC 
state Eq. (ITT1) to the one proposed in Ref. , which comes from the computation of the amount of gluon production, 
starting from a dense-dilute configuration of initial protons. 

Let us apply our approach to Gaussian models of unintegrated gluon distributions, in which we introduce a pa- 
rameter generalization of the GBW model [§] in order to clarify the contributions to the entropy. The GBW model 
indeed is used by Q for its definition of entropy which we wish to compare with our formulation. We thus consider 
unintegrated gluon distributions ip(k,Y) ([1} with a Gaussian tail, namely 

V{k,Y) d 2 k^P(u = R 2 k 2 ) du = 4>(u) e~ u du 

~ r -1 (/e) u"- 1 e~ v du , (12) 

where the non-Gaussian prefactor i/)(u), describing the behavior in the low- A: range of gluon momentum, is more 
simply parametrized by a power-like term u 1 ^" 1 . We shall call k the "overlap parameter" for reasons becoming clear 
later on in this section. Note that k = 2, which is the GBW value, corresponds to a dipole distribution verifying the 
transparency property at zero dipole size 

Definition (|12[) is indeed a slight generalization of the unintegrated gluon distribution coming from the Golec-Biernat 
Wusthoff dipole model [7|, |9| obtained for the value k = 2. 

Starting by implementing (|12l) in the dynamical entropy definition ([5]), one gets 

and, more explicitly for the K-dependent family (1121) 



iYi->y2 _ 



k 2 (R\- Rl) - k log 



R 2 



V{k 2 R\) R\dk 2 



= (k 2 ) Y2 [Rl-RD-K log (R 2 /R 2 ) = k { (Rl/R 2 2 - 1) - log (Rj/R 2 ) } > . (14) 

One can easily check the strict final inequality (fT4")l for all Y2 > Y\, i.e. R2 < R\ while it goes smoothly to zero (as 

k/2 (R\/Rl - if) when i? 2 ->i?i- 

Let us list relevant properties of formula (|14[) . 



6 



• Dynamical entropy of the CGC state. Using the relation dill) , we can compute the total entropy density 
dS ^ 2 /nR^dy for a CGC state of overall transverse size Rt- Assuming the framework of an initial distribution 
of independent transverse "cells" with size R\ , it is given by 



1 dS Y ^ Y - 



1 1 



nR T 



dy 



Air 2 N r 



R\ 



N 2 



1 



4ir 2 N r a< 



Qi - Qi 1 + lo 



Q 



(15) 



expressed in terms of the saturation momentum dependence Q s = 1/R S . Note that at large Q^/Qu formula 
(1151) shows that the target dynamical entropy by unit of transverse area and rapidity becomes independent of 
the initial condition at Y\ and is directly proportional to Q 2 . The logarithmic corrections provide a smooth 
transition to quasi reversibility dS^-0 when Q\jQ\-^\. 

• "Microscopic" vs. "macroscopic" entropy: determination of thea parameter. Let us now derive the dynamical 
CGC entropy density precisely for the case considered in Ref.[7|, and check the matching with the expression 
obtained in [7| from gluon production in the dense-dilute collision case. We tend to call it "macroscopic" entropy 
since it relies on a derivation from thermodynamic concepts. In our case, that we call "microscopic" we proceed 
from a counting of the configurations in transverse momentum phase-space parametrized by the probability 
V{k,Y), from (|T2]). 

For p—p scattering, considered in 0, we use, for simplicity, ([T5j) for the GBW gluon distribution with k = 2. We 
find for the dynamical ("microscopic") entropy 



1 

A T 



dS _ 2fi(N 2 - 1) 
dy Tr 2 a s N c 



1 + 



Vl 

Q 2 Qi 



(16) 



with k — 2 and the proton mean transverse area At = itR T /A. 

This results correctly matches the one 3 obtained in the derivation of ref. 0] which also gives the same behavior 
as the first term in formula (|16[) . independent of the initial condition. Hence the two derivations of the entropy 
of a CGC medium are consistent with each other despite their marked difference. Indeed, one comes from the 
gluon production process in a dilute-dense collision and the other is based on the configurations of the CGC 
medium itself. 

Moreover, for an exact matching between the two formulas, one is led to identify 



3tt 
T 



(18) 



which gives a non trivial estimate of the number of effectively independent gluonic degrees of freedom within 
a given color correlation length R s . Interestingly, it comes as a pure number, which may be characteristic of 
the average number of gluons inside a colorless cluster of a CGC medium generated by saturation. It would be 
worth checking that this remains valid for more elaborated UGD's. 

• Strong irreversibility of the formation of a CGC state. It is interesting to note that, using the same method, one 
can compute from © the dynamical entropy Y, R2 ~^ Rl associated with the backward process , i.e the expansion 
R2 — >Ri > i?2 ■ By exchanging R 2 -n- R\ in (fT4|) one finds 



£,R2->-fli 



lot 



P(k,Y! 
V(k,Y 2 



= n{\og(R\/Rl)-{l-Rl/R\)} >0 



(19) 



It is easy to check also in that case the strict inequality for R\ > R 2 , with a same smooth behavior towards 
the equality as R 2 ^-Ri. However, one finds a logarithmic behavior log (Ri/R 2 ) by contrast with the quadratic 
power- like behavior (|14p . In fact, in that case the main contribution to the dissipative work is due to the free 
energy term in (|10[) . while the "work" contribution is almost zero. The interpretation is that (c/. [H[) only 
few gluons within the stochastic distribution of transverse momenta feel the expansion of the color correlation 



3 The relevant formula (25) in [t| reads, in our notations 

1 dS 3(7V C 2 -1) 2 

-Q B +cst. , (17) 



At dy irasN c 



7 



length. Hence the "gluon compression" mechanism is in general much more dissipative-and thus irreversible- 
than "gluon expansion" . This is even more remarkable, when assuming a "cyclic" process, namely a compression 
R% — > i?2 followed by the reverse expansion R\ — > R2. The overall balance of entropy AE^(i?i, R2), following 
Eqs. (|14ll9l) writes 

A^(R 1 ,R 2 ) = K [ \ 2 £ ) (20) 

where the last approximation is for Ri ^> R 2 large. Formula (f2T)|) makes explicit the irreversible feature of the 
color correlation size changing processes 4 . 

From this, one can induce that the QCD saturation mechanism due to the high density of gluons is generating a 
large entropy and thus highly dissipative in thermodynamic language. Hence the dynamical entropy ([5]) seems to 
capture interesting irreversible physical features of the formation of dense QCD states, and eventually, through 
the cyclic process, their relaxation towards equilibrium. 

• Comparison with a 1- dimensional ideal gas model moved by a piston. It is quite instructive to compare the 
results for the QCD two-dimensional dynamical entropy ((HJ) with the non-relativistic one-dimensional ideal gas 
model of Ref.[18]. One considers the rapid action Li t» L / < Li of a piston on a longitudinal box containing an 
ideal gas of independent particles, starting with temperature T. One can compute the work distribution due to 
the piston either by expansion [18| or compression [24fl The piston model is a well-known far-from-equilibrium 
example where the Jarzynski identity (j9|) can be explicitly verified, and one finds for the average dissipative 
work (W) — AF (and thus the dynamical entropy) 

Af ^l 2 1 K Li/La _ 1 _ i og (L 1 /L 2 ) 

AF L2 ^ Ll } AF L2 ^ Ll oclog(Li/L 2 ) ■ (21) 

It is tempting to interpret our results (fT4|) and (fT9j) as the radial analogue of the two equations of (|2"Tj) . In some 
sense, the QCD field-theoretical saturation mechanism can be thought of as due to multiple shocks of individual 
gluons on the effective "walls" due to the shrinkage of the color correlation length, as a mean-field effect on an 
individual gluon reproducing the many-body interactions of gluons at saturation. Rapid expansion is much less 
entropy producing, since it creates almost no work, leaving only, as a dissipative mechanism, the logarithmic 
free-energy cost — AF. As a side remark, it is interesting to figure out that saturated gluons within a small color 
correlation "cell" may behave as non relativistic particles, if they acquire an effective mass M ~ Q s as suggested 
in0. 

• Role of the "overlap factor" k. The factor n, which parameterizes the low-fc range of the gluon transverse 
momentum distribution (cf. (|12p ). appears as an overall multiplicative factor in our resulting formulas for the 
dynamical entropy (|14I15I19[) without modifying the quadratic dependence in Q 2 . It is linked to the structure 
of the gluon probability distributions V(k,Y2) vs. V(k,Y\), cf. (jl"4")) . Indeed, a smaller (resp. larger) value of k 
corresponds to wider (resp. less) overlap between the two distributions and thus lower (resp. higher) amount of 
dynamical entropy. It appears to be much dependent on the low-p^ range of the UGD's. 

A remark here is adequate concerning the Weiszsacker- Williams gluon distribution, which is also discussed in the 
framework of UGD's (25J. The Weiszsacker- Williams gluon distribution has the physical interpretation as the 
number density of gluons inside the hadron/nucleus in light-cone gauge, but is not appearing in the cross-section 
observables (25[. The main difference with the one used in (TTSj) (and also in 0]) is the infra-red logarithmic 
behavior for k 2 —>0, would correspond to an overlap factor n s» 1 in (fl"4"|) . In our definition (fT5|) of the dynamical 
entropy, we choose the "dipole gluon distribution" and formula (|15l) for the following reasons: it corresponds 
to a suppression of the color radiation at null transverse momentum, and a neat maximum at the saturation 
momentum Q s , which is in a good analogy with the classical notion of a gas of particles in a "box" of size 
R s = 1/Qsi which makes consistent the aspects of classical far-from equilibrium statistical physics on which 
relies our definition. This is related to the fact that QCD dipoles are colorless states which follow classical 
dynamics, at least at leading logs order of the perturbative expansion. 



s Ll ^ L2 = i{(w>- 



4 Such "cyclic" process could be eventually describing the entropy production obtained through the formation and subsequent relaxation 
of a CGC state, when evaluated at weak QCD coupling. It would be interesting to relate it to the entropy production through 
thermalization. 
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IV. EXTENSION TO THE DENSE INITIAL STATE OF HEAVY-ION REACTIONS. 



The pre-equilibrium state of matter reached immediately after the heavy-ion collision, is only a transient one. In the 
most common scenario, it emanates from the interaction of the two CGC states representing the wave-function of the 
incident colliding particles, forming "at time zero" what is thought to be in the literature the glasma phase Q. Then 
it evolves towards the hydrodynamic QGP phase, during a strikingly short-time period, as indicated by simulations, 
which is referred to as plasma thermalization. This transition is still the less understood part of the QCD plasma 
formation in heavy-ion reactions. We shall investigate how the notion of a dynamical entropy may provide some new 
light on the questions raised by the rapid thermalization found in heavy-ion collisions. 

The two incident nuclei considered in the center-of-mass frame are boosted and thus can be described as two 
independent CGC states. Hence, our formalism may in principle be applied to the initial CGC states of both nuclei. 
However, one has to take into account the glasma state formed from these CGC states at the interaction time. 
Hence, the application of our method defined in the previous section requires some extension. In particular, the 
significant matching with the "macroscopic entropy" of Ref.Q, see previous section, was made considering "dense- 
dilute" collision while the glasma state corresponds to a "dense-dense" configuration in the overall center-of-mass 
frame. In the following, we shall make a proposal how one can extend the notion of dynamical entropy for the glasma 
and its impact on the thermalization problem. 



• Dynamical entropy of the glasma. 

Let us consider the glasma phase of Ref. |28l| , as an application example of our approach. Extending the dynamical 
entropy concept to the glasma defined in [28j] appears nontrivial, since one has to consider a classical Yang-Mills field 
calculation taking as initial value the result of the JIMWLK evolution for each nucleus. Those calculations lead to a 
prediction for the gluon distribution in the glasma 4> g iasma , y) where p± is the transverse gluon momentum and 



y, t he rapidity evolution away from the initial condition at y = , taken in [28| as the McLerran-Venugopalan model 
23]. What is numerically observed is that the resulting glasma spectrum satisfies geometric scaling using the gluon 
(i.e. QCD adjoint) saturation momentum Q S (Y). We are thus led to apply a similar recipe as in the CGC case, with 
the difference that new nonlinear effects due to the classical Yang-Mills calculations are present, especially in the low 
p± range. One thus may proceed as follows. 

Defining a probability distribution, similarly to ([3]), 

V glasma (p ± ,Y) d 2 P± = ^srnaip^Y) dPk ^ p = 2 2 

J P glasma (P±,Y) d 2 p± 

one is led to propose for the glasma dynamical entropy density (in a transverse cell of order R ~ 1/Q S (Y), for the 
evolution Yi —> Y 2 



^Ml°g p7'T ' y > 0. (23) 

\ / glasma\P±, II) / y 2 

Assuming for simplicity a Gaussian from of the distribution, we consider a rapidity evolution of the glasma state 
in the range where geometric scaling with saturation momentum Q S (Y) is valid. Taking also into account in a 
phenomenological way the observed depletion [28j of gluon radiation in the small range of p±, one may write a 
formula similar to (fT2")l being used for the CGC states in section Hill namely 

V glasraa ^ = ^ {p±/Qs) e - P l/Ql # p± „ _L_ ^ ' e ~A/Ql ^ p± ^ 

where the last approximation comes from considering the same input as in Ref.Q, i.e. (I12p and the small- 
p_i_regularization prefactor ip(px/Qs)) being approximated using an appropriate overlap factor K g i. Applying now 
our definition ([5]) we get 

^laZla ~ HI QllQl > (25) 

where, as usual Qi = Q s (Yi). It would be interesting to make a more precise evaluation of the glasma dynamical 
entropy and compare it with the sum of entropies of the initial CGC states for the same rapidity increase. This 
require to go beyond the Gaussian toy model we are using and we postpone this to further study. 



9 




FIG. 1: Piston model: Compression and Thermalization. Left: Thermal state at temperature To and length Middle: Far- 
from equilibrium compressed state. Right: Thermalization through entropy production to the heat bath and relaxation to a 
denser thermal state at temperature T and length Lf. 



• Entropy production and thermalization in heavy-ion collisions 

The pre-equilibrium state of matter reached immediately after the heavy-ion collision, is only a transient one. It 
evolves, during a quite short time as indicated by simulations, towards the hydrodynamic QGP phase which will give 
rise finally to hadronization. This stage, driving the interacting medium from the CGC initial nucleus states through 
the glasma formation to the hydrodynamic behavior of the QGP, called thermalization is not well understood, due to 
lack of a microscopic theory of thermalization based on QCD. It may require a strong coupling treatment of collisions 
in QCD which is not yet theoretically available. Some hints have been obtained [|| from the AdS/CFT correspondence 
for the similar problem formulated in the case of the maximally symmetric J\f — 4 gauge field theory. 

In this context, let us examine what can be obtained from the formulation of the dynamical entropy for the 
initial state before thermalization, using the previous results. For this sake, we will confront both macroscopic and 
microscopic points of view. 

From a macroscopic point-of-view inspired by out-of-equilibrium statistical physics, there exists an intimate relation 
between the dynamical entropy of an initial state as defined in our approach and the entropy production during 
thermalization. This has already been discussed in the case of the dilute-dense configuration of the initial states, see 
section Hill Or goal now is to extend it to the dense-dense case which is typical for heavy- ion collisions in the central 
rapidity region. Second we will consider a microscopic approach, namely the strong coupling scenario of Ref.Q, 
where the notion of an "initial entropy" has been introduced using a dual gravitational version of the Af = 4 gauge 
field theory. We will conjecture that its determination for QCD may be based on the dynamical entropy of the glasma 
state. 

Coming back for this sake to the Jarzynski and Hatano-Sasa identities. A key feature is their interpretation in 
terms of the entropy generated by a fast non-equilibrium process in a finite time, and released to the environment by 
relaxation. In the case of the Jarzynski identity §§§ , see Appendix [B] and assuming a thermodynamic coupling with 
an outside heat bath at temperature T, the exponent is related to the amount of average dissipative work (W dt ^), 
itself connected to the production of entropy AS* by the relation 

AS=±(-W diff ) = ±((W)-AF)>0. (26) 

In practice, the produced entropy comes from the relaxation of the out-of-equilibrium state by freezing the external 
force acting on the system, which finally towards thermal equilibrium at the initial value of the temperature T. 

This process is, among others, well exemplified by the exact results obtained for one-dimensional piston model [l8j . 
see Fig.l. Starting from the thermal state of an ideal gas at temperature T in a box of initial length Li, the action 
of the piston, with speed u during a short time Ati, shrinks the box size till a final size Li~^Lf=L— |w|Ati which 
is then kept fixed, leaving the system relax towards the initial temperature 5 . Indeed, the action of the piston gives 
rise to the far-from-equilibrium state of Fig.l middle, with the particles of gas being far from thermal equilibrium 
since most part of them acquires momentum and kinetic energy by shocks on the piston [l8j . Then its thermalization 
through the production of entropy makes the system relaxing to the thermal state at temperature T and length Lf, 
depicted in Fig.l right. 



Another interesting exact model is the adiabatic expansion and compression of a dilute interacting gas |17|I 
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It is interesting to note that the estimate (|25|) of the dynamical entropy of the initial glasma state appears as a 
simple extension of relation (|26[) with 

AS = I (w*//) ^ e«^« > ; i(W)^/%z 5 i - log ) , (27) 

for a two-dimensional version of the piston model. The physical interpretation is that the evolution of the glasma 
state is obtained from the compression of gluons confined within a box of color correlation length Ri — > R 2 due to 
the increasing rate of branching and recombination of gluons at higher density. 

The precise determination of AS asks for a mechanism incorporating the expansion of the QCD matter, the 
thermalization process, the subsequent hydrodynamic evolution and eventually the hadronization of the QGP. While 
such a theoretical treatment is still unknown for QCD (beyond phenomenological approaches), a refined treatment 
of thermalization through the hydrodynamic evolution could be obtained in the strong coupling scheme of Ref . |3| for 
the M = 4 SYM theory. Let us then consider this case in the context of our approach. It has been shown that the 
boost-invariant proper-time expansion of the Af =4 SYM plasma leads to a hydrodynamic stage [3] with small viscosity 
over entropy ratio. Numerical solutions of the early-time evolution using the gravity dual allows one to examine the 
properties of the pre-equilibrium evolution, starting with various initial conditions which are not specified in [U , and 
thus not related to the glasma states. We try and fill this gap in the sequel. 

The main common feature of the solutions of the strong coupling studies using the gauge/gravity duality, beyond 
the variety of initial conditions, are the short proper-time evolution towards the hydrodynamic stage, before which 
they differ. Those regularities, such as entropy variation and the thermalization temperature seem to only depend on 
the initial energy and entropy densities, namely 



e(0) = N*- ^.^(0) ; ^ =iV c 2 .i 7 r 2 .T e 2 // (0). S (0) (*» 



where one introduced an "initial effective temperature" T e ff(0) and an "initial entropy" density s(0) by degree of 
freedom, both defined by (|28p from the corresponding energy and entropy densities by unit of space-time rapidity r\ 
and transverse area At- 

Note that the initial state is not at equilibrium and thus T e //(0) is not a physical temperature and s(0) is not the 
usual entropy at equilibrium. Following a notion of "dynamical entropy" in general relativity, it is defined via the 
holographic correspondence as the element of the apparent horizon's area situated at geodesic distance at proper-time 
t = from the origin in the bulk space (see Q and references therein). 

Factorizing the color degrees of freedom, we note that (|28|) gives rise to the following relation 

1 (3/4) dS/dn , n , 



8(0) ■ T-/J0) , (29) 



A T e (0) ~ Jyu > 

valid for each degree of freedom of the initial state. Note that the correction factor 3/4 is coming from the known 
ratio (see e.g. 0, section 3.6) between the strong coupling entropy at finite temperature and the ideal gas one. We 
shall not take into account this factor in the following, which deserves more study going beyond our context. 

Considering now the corresponding densities for dense QCD initial states, and for application, the Gaussian model 
(f2"2"|) . we may write formulas in parallel with (|28l) and (|2"9"|) . For this sake, one starts from the formulas for the glasma 
energy density 6 and for the dynamical entropy (|25j) . with the standard counting of gluon degrees of freedom (assuming 
for sake of simplicity the same occupation number l/(4ir 2 a s N c ) in rapidity as for CGC states, cancelled in the ratio). 
One obtains 

^c-l r* 1 dSgl N?-l -2 

€sl = ^^N c - Q v i 5 ^^ = K9 '-^-47wv;- Q ^' (30) 

where e g i is the initial 7 energy density of the glasma, 1/Q g i = R g i is the size of the glasma color correlation length, 
K g i the overlap factor in the Gaussian case (I24[) and (i g i is the number of gluonic degrees of freedom in the color 



6 For simplicity we use for the energy density in rapidity the standard factor. For a more refined determination, see |27||. 

7 There is a subtlety in the determination of e g ; if using as an input the Weiszsacker- Williams definition of the UGD [27| which diverges 
at t = 0. We here still use for consistency of the Gaussian approach the "dipole" UGD. 
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correlation "cell" of the glasma (cf. (|18|l for its determination for the CGC). One finally note that (|30|) gives rise to 
the following normalization independent relation 

1 dS g i/dy 2 

- Kgl ■ n gi ■ Qni ■ ( 31 ) 



At e g i(0) 



Matching (|3Tj) with ([29]) can be performed, since the proportionality of T e ff(0)/Q g i is expected from the scale 
invariance of both the N = 4 SYM theory (from conformal invariance) and the dense state distributions in QCD at 
saturation (from geometric scaling). 

Now, in order to get an estimate of the initial entropy s(0) taking as an example the CGC case 0, one uses K g i = 2, 
with T e ff(0)/Q g i — l/(27r), from the determination of the temperature [29| through the Unruh property. In the 
Bjorken boost-invariant regime ([l|, first reference) one can identify the space-like and energy- momentum definitions 
of rapidity r) = y. Finally we assume the same value \i g i = 3ir/2 as in (|18[) for the number of gluons by color cell. We 
then find 

s(0) - n gl ■ (x gl • (T eff (0)/Q gl ) 2 => J- . (32) 

where the last number comes from the Gaussian form of a CGC state, taken as an example of the dense QCD state. 
Note that this number is in the ballpark of the range < s(0) < 1/2 obtained in [3J. It corresponds to a quite sizable 
entropy density. 

More generally, the key order parameter s(0) of the thermalization analysis in J\f = 4 SYM theory is in direct relation 
with the "overlap parameter" K g i, characterizing the strength of the dynamical entropy of the glasma state (in TV" =4 
SYM). This is consistent with the entropy increasing with decreasing overlap of probability distributions. 

Recalling the main results of the numerical analyses showing Q that the system evolves quite rapidly towards a 
hydro dynamic behavior, even before isotropization occurs, the rate of entropy production Sf— s(0) and thermalization 
proper-time (indeed, better to tell "hydrodynamization time") T t h in units of T~^(0) are driven by the value of s(0). 
In our corresponding QCD evaluation (|32[) for the Gaussian model, it is the value of the overlap parameter k which 
is relevant. Higher is the dynamic entropy parameter k, i.e. smaller is the overlap between the gluon probability 
distributions at Q2 vs. Q%, where R± = 1/Qi is taken for convenience of order of the standard color correlation 
length at rest. Then, greater is the initial entropy s(0) and the subsequent increase of entropy Sf— s(0) during 
thermalization. Following [3[, the thermalization time is consequently shorter. Hence, in our framework, it means 
that the thermalization time depends on the "overlap parameter" , and thus on the strength of the dynamical entropy 
of the initial dense saturation state. 



V. CONCLUSIONS 



Let us give our conclusions. We have proposed a definition of a dynamical entropy for dense QCD states of matter at 
high energies. This definition is based on statistical physics tools which have been shown to be valid for large classes of 
non-equilibrium processes. Our approach is inspired by the Jarzynski and Hatano-Sasa fluctuation theorems [HI, [l6| 
and is applied to the formation of dense QCD matter through an energy-dependent evolution scheme defined in the 
QCD field theoretical framework at weak coupling. It describes the evolution of the color correlation length Ri — >B,2, 
corresponding to the rapidity evolution Yj — > Y2 of a dense QCD medium of gluons which is the result of gluon 
branching and recombination. 

The derivation of an identity (O provides a guide for the definition of a dynamical entropy functional verifying 
positivity, similar to the 2 nd principle of irreversible thermodynamics. It acquires an interesting physical picture 
when the gluon distribution of the dense QCD medium reaches the geometric scaling regime. In that stage, the 
boost-invariant distribution of gluons can be described as transverse cells whose radial size is constrained by the 
collective effect of the others to be the saturation scale R s = 1/Q S (Y). The dynamical entropy is expressed as an 
overlap functional between the gluon distributions at different total rapidities Y\ and Y2 > Y\ and thus saturation 
sizes R„(Yi) and R S (Y 2 ) < i? s (Yi). 
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Summary 



Summarizing concretely our results: 
• Dynamical Entropy of CGC states 



It is possible to define an entropy functional (by unit gluonic degree of freedom) 

where V{k, Y) is the probability distribution of gluon transverse momentum at rapidity Y. It is defined using the 
QCD dipole unintegrated gluon distribution after normalization to unity. £ Yl_>Y2 is the dynamical entropy density 
(by transverse space "cell" ) during a far-from-equilibrium "compression" process during which the correlation length 
Ri shrinks to R2 < R\. 

In the geometric scaling set-up, the color correlation length of transverse size R s appears explicitly as 

SMl = ('° g §¥j). - be 3~f {W—AF— } . (34) 



Hence, (|33[) in the geometric scaling form (|34[) admits the interpretation of a genuine entropy in a thermodynamical 
framework, see Appendix [XJ since the dissipative work performed during this process corresponds to the difference 
between the average work (\ogV{R\k 2 )/T'(R\k 2 yj Y and the variation of the free energy \ogR\ fR\ between the initial 
and final "cells" , during the "compression" process R\ — > R2 due to the increase of the gluon density when the total 
rapidity evolves Y\ — > Y2 . 

In a non-thermodynamic framework, see Appendix [B] where the CGC states are assimilated to non-equilibrium 
"stationary states" , then Y, Rl ~^ R2 plays the role of a dynamical entropy acquired by the CGC "initial stationary state" 
from a rapidity evolution and released through relaxation. 

• Application : the Gaussian model. 

As an example, we considered dipole gluon densities of the form </>(fc,i?i) ~ (kRi) 2 ^ -1 ^ exp (— k 2 R 2 ^j generalizing 
the GBW model @, for which one explicitly finds 

S fl ^=«{(Ql/Q?-l)-log(Ql/g?)} >0, (35) 

where Q s (ii,l2) = 1/J?n ; 2)) *- e - the saturation transverse momenta at the initial and final stage of the rapidity 
evolution. The parameter re is an overlap factor characterizing the amount of overlap between the initial and final 
distributions, especially at low transverse momentum. Greater is re, smaller is the overlap, and thus stronger is 
the dynamical entropy density generated by the evolution Yy — »-l2 > due to the increasing number of branching and 
recombination QCD vertices involved in the process. 

Interestingly, the result ([33)1 was proved to be identical for large AY, up to a factor fi = 3tt/2, to the evaluation 
of the entropy variation of a CGC state coming from a macroscopic formalism, i.e. a thermodynamic interpretation 
Q of the distribution of gluons produced in a dense-dilute collision. The characteristic factor /i can be interpreted 
as the (average) number of gluonic degrees of freedom inside a transverse CGC "cell". To our knowledge, this is a 
new result coming from the comparison of two different but complementary approaches to the notion of entropy for 
a CGC medium. This appears to be in conformity with the statistical physics interpretation, via the comparison 
between microscopic vs. macroscopic approaches, which relates the work fluctuations -here given the momentum 
distribution in the CGC state- to the entropy production - here, the gluon multiplicity produced in a "dense-dilute" 
collision probing the structure of the dense medium. 

• Dynamical entropy of the glasma in heavy-ion collisions. 

We considered the class of models using the Color Glass Condensate formalism describing the initial nuclei at the 
collision time and leading to the glasma as the initial interacting state of heavy-ion collisions. Using the rapidity 
evolution of the glasma distribution of gluons, as e.g. performed in Refs. (27l. l28j . it is possible, at least approximatively, 
to define the dynamical entropy for an rapidity variation 8y in a way similar to the CGC states. 

As an illustration, we used a phenomenological Gaussian model for the probability distribution, namely 

where, from the analysis of (28j, one verifies geometric scaling with momentum scale Q 2 corresponding to the adjoint 
QCD representation, and a regularization at smaller momentum due to the nonlinear effect of the classical Yang-Mills 
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fields initiated by the initial CGC states. It is described here by a phenomenological value n g \ of the overlap factor. 
One then gets in this example at leading order, for one gluonic degree of freedom, 

• Hints for thermalization: 

The evaluation through QCD calculations of the dynamical entropy for the initial state of heavy-ion collisions allows 
for interesting hints on the problem of the transition towards the thermalized hydrodynamic stage of the Quark-Gluon 
Plasma. 

Noting that the expression (|37| for the dynamical entropy density of a glasma state can be a candidate for the 
"initial entropy" density s(0) appearing as a key parameter in the strong coupling thermalization mechanism Q of 
the ADS/CFT correspondence, one finds the relation 

s(0)-^-AV-{^^} 2 , (38) 

where T e ff(0)/Q g i is the ratio between the effective temperature (determined by the initial energy density) and the 
saturation scale of the glasma, expected to be a universal constant, independent of the total rapidity Y. The parameter 
H g i is the number of gluonic degrees of freedom in a glasma "cell" . Assuming for simplicity the same values for k, fi 
as for the Gaussian CGC distribution, we find values of s(0) in the range compatible with the AdS/CFT framework 
ofRef.0]. 

It is interesting to see that (|38j) gives a realization of the main physical idea behind the Jarzynski and Hatano-Sasa 
identities, namely the connection between the work (and thus momentum) probability distribution and the overall 
produced entropy by a far-from equilibrium process. Indeed, a large (resp. small) overlap factor K g i means that 
the gluon momentum variations induced by the rapidity (and thus energy) evolution of a glasma state have been 
large (resp. small). Hence the contribution to the subsequent entropy (and thus particle) production is larger (resp. 
smaller). 

In practice, and if the input parameter values are taken from the CGC example in QCD and the AdS/CFT example 
from in Af = 4 SYM theory, we find a rather large initial entropy density, in the range considered in [3|] . It remains 
to see whether more realistic applications, both for QCD and AdS/CFT confirm this first approximation. 



Outlook 



Many pending questions about non-equilibrium processes occurring during a high-energy heavy-ion collisions which 
may be discussed using the tools of non-equilibrium statistical mechanics, as we started to perform in the present 
study. Let us mention some of the prospects in that direction. 

• Dynamical Entropy of dense QCD states : theory 

Considering our application to QCD far-from-equilibrium processes, the property of the stochastic identities (e.g. 
[HI is the possibility of finding an expression for the entropy -usually related with the final multiplicity- 
using the knowledge of the momentum distribution in the initial state. Hence, on the theoretical level, the question 
arises how the properties of the dynamical entropy are related to the basics of the QCD evolution equations in the 
saturation setting, namely the BK or JIMWLK equations, see e.g. Q. Indeed, they are known to lead to geometric 
scaling gluon distributions when the rapidity Y increases, which is used in our derivation. One more specific question 
is the role of the geometric scaling violations, which are known to occur in the full solution of the QCD equations 
[n| . Hence, one may ask what is the mathematical meaning of the general expression of the dynamical entropy 

Ss / log gy^\ >„, (39, 



Vlk.Y,) 

where V(k; Y^) is proportional (with a suitable normalization) to the derivative of the gluon structure function F(k, Y), 
solution of the BK/ JIMWLK QCD evolution equation as a function of the rapidity Y. This can be addressed as a 
question about integro-differential equations of F-KPP type which admit asymptotic traveling wave solutions, and 
are the proper mathematical formulation of geometric scaling [ll[ ■ 
• Dynamical Entropy of dense QCD states : phenomenology 

On the phenomenological basis, it would be interesting to make a more quantitative estimate of the various entropy 
contributions during all expected successive steps of a heavy-ion reaction using the phenomenological knowledge one 
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has acquired about these states. Namely, we can usefully consider physical parameterizations of CGC states and of 
the glasma phase, and of their evolution towards the hydrodynamic phase. These problems deserve future studies. 

Finally, we may hope that the gap of knowledge yet existing for a proper general and operative formulation of 
nonequilibrium processes in quantum field theories may be reduced thanks to studies on heavy-ion collisions. We 
conjecture that the existence of generalized work identities may give some hints for elaborating such a general 
framework, may be in relation with the formalism of Wilson lines and loops, which is rooted in the bases of QCD, 
both at weak and strong coupling. 

• The thermalization problem 

The formalism and calculation presented in section IIVI for the initial state entropy of an heavy-ion collision are 
in fact the weak coupling analogous of the strong coupling ones in the framework of the AdS/CFT correspondence 
performed in Ref.[3J. A comparison between both approaches, see e.g. (|38p . shows that a matching can be obtained 
between the glasma initial state 8 and the subsequent evolution within the strong coupling scheme, the question now 
arises to eventually make a more direct comparison between the definition of entropy in those dynamical regimes: the 
one coming from the Shannon-type definition (|33[) and the one coming from the dual gravitational theory related to 
the local area of the apparent horizon [3j . Our formalism may also be useful in the thermalization approach based on 
instabilities, e.g. Ref.[30|. 



Appendix A: The Jarzynski identity 

Generic identities have been derived when the evolution is fast and thus very far from equilibrium. These are the 
Jarzynski [13J and Hatano-Sasa [16( (also Crooks [15j) fluctuation relations, which we now briefly describe in the text. 

It is well-known from statistical physics textbooks that inequalities involving the average work (W) performed by 
a statistical system verifies thermodynamic inequalities related to the second principle of thermodynamics, namely 
the increase of the total entropy for an isolated system. For instance, for the thermodynamic transition between two 
equilibrium states, one automatically verifies the celebrated second principle of thermodynamics 

AS = i«W) - AF) > = i x (W Diff ) > , (Al) 

where AS is the increase of entropy, T, the temperature, and W is the work performed in the medium i.e. AE = (W) 
the total energy change, AF the variation of the free energy of the system. In such a relation, W — AF = Woiss 
is the dissipative work paid by the system during a nonequilibrium process. The equality to zero in (|A1[) is only for 
reversible processes, through slow dissipation-less evolution. However, as in particular shown in by the work identities 
[HI there are necessarily fluctuations with Woiss < in the dissipative work distribution [l3| . 

To be a bit more precise in this short introduction, the Jarzinsky identity relates the stochastic distribution of 
thermodynamical works in the process A — > B to the free energy balance AF between the two equilibrium states 
A— > C. Interestingly, the amount of dissipative work y\? Dl ff = W — AF during the process A — > B is then related to 
the entropy production AS = 1/T((W) — AF) > , if the state B is able to relax towards the temperature T, keeping 
the driving parameter constant. 

The Jarzynski identity [l]| reads 

(e- w l T ) = e- AF / T . (A2) 

It is to realize that, while verifying the second principle (W) — AF = AS > thanks to the Jensen identity f2ll . 
it implies deep properties on far-from-equilibrium fluctuations. Indeed, it can be shown in many examples (cf. |14| |) 
that a dominant contribution to (|A2[) is given by rare out-of-equilibrium fluctuations very far from those dominating 
the average (W) and thus AS* in ([AT]) . 



Appendix B: The Hatano-Sasa identity 

The nonequilibrium identities have a much larger range of applicability than the one yet considered initially [l3j . 
Among the various extensions, we will select for our purposes the one |16| which applies for purely non-equilibrium 



We expect the glasma state to be similar QCD and Af = 4 SYM theories, which are both scale-invariant in the leading-log perturbative 
domain. 



15 



systems without reference to the notion of temperature. For this sake, one considers a situation where one has a 
system governed by a dynamical parameter A, such that at each value of A, it corresponds to a stationary phase 
space spectrum with probability distribution V at (z; \^ z ^ w here the variable z describes the phase-space. Then the 
following equality holds 16] 



exp - 



, dX d\ogV(z;X,T) 
dT Tr dX 



— dz exp - 



T2 , dX d\ogV{z;X,T) 
, dT Tr dX 



P(z;X 2 ,t 2 ) = 1 



(Bl) 



where by definition V(z;X, r) = V at (z(T);A(r)) is the stationary solution for the value A(r) in the "frozen" phase 
space variables at time r. Note that the resulting identity is a -priori not dependent on the arbitrary "history" 4^ of 
the non equilibrium mechanism. 

In fact, the quantity — log'P Stat (2;; A) provides a candidate definition for a quantity playing the role, after integration 
in phase-space of a "dynamical entropy " denoted S. Indeed, for the dynamical non-equilibrium process, where one 
starts with a state corresponding to Ai to a state A2, the increase of 'dynamical entropy" gives rise to the following 
formula 



AS 



Ai->A 2 _ 



= log 



^ btat (Ai) 



= - / log 



jStat 



(z, Ai 



V Stat (z,X 2 ) 



x 7> stat (z;A 2 ) dz 



(B2) 



In fact, our approach, see (O, Eq. dBip is simply obtained by noting that there is no parameter dependence on time, 
but only via X Y which allows for an explicit integration in (|B2[) . Moreover, geometric scaling ensures the identity 
(jBljl since "P stat (z(T); A(r)) depends only one one variable u(t) = z(X(t)). 

As can be noticed, a key feature is that no equivalent of a temperature appears in the relation (|B1[) . In the Jarzynski 
case using Boltzmann factors, where temperature appears explicitly, the Hatano-Sasa identity is proven to reduce to 
the original Jarzynski equation (|A2|) . namely 



e -^(W-AF)\ 



1. 



(B3) 



where W dlss is the dissipative part of the work. 
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